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ABSTRACT 

Context. For the past two decades, the three-dimensional structure of sunspots has been studied extensively. A recent improvement in 
the Stokes inversion technique prompts us to revisit the depth-dependent properties of sunspots. 

Aims. In the present work, we aim to investigate the global depth-dependent thermal, velocity, and magnetic properties of a sunspot, 
as well as the interconnection between various local properties. 

Methods. We analysed high-quality Stokes profiles of the disk-centred, regular, leading sunspot of NOAA AR 10933, acquired by the 
Solar Optical Telescope/Spectropolarimeter (SOT/SP) on board the Hinode spacecraft. To obtain depth-dependent stratification of the 
physical parameters, we used the recently developed, spatially coupled version of the SPINOR inversion code. 

Results. First, we study the azimuthally averaged physical parameters of the sunspot. We find that the vertical temperature gradient 
in the lower- to mid-photosphere is at its weakest in the umbra, while it is considerably stronger in the penumbra, and stronger still 
in the spot’s surroundings. The azimuthally averaged field becomes more horizontal with radial distance from the centre of the spot, 
but more vertical with height. At continuum optical depth unity, the line-of-sight velocity shows an average upflow of ~ 300 ms"^ 
in the inner penumbra and an average downflow of ~ 1300 ms“^ in the outer penumbra. The downflow continues outside the visible 
penumbral boundary. The sunspot shows, at most, a moderate negative twist of < 5° at log(T) = 0, which increases with height. The 
sunspot umbra and the spines of the penumbra show considerable similarity with regard to their physical properties, albeit with some 
quantitative differences (weaker, somewhat more horizontal fields in spines, commensurate with their location being further away 
from the sunspot’s core). The temperature shows a general anti-correlation with the held strength, with the exception of the heads of 
penumbral filaments, where a weak positive correlation is found. The dependence of the physical parameters on each other over the 
full sunspot shows a qualitative similarity to that of a standard penumbral filament and its surrounding spines. 

Conclusions. The large-scale variation in the physical parameters of a sunspot at various optical depths is presented. Our results 
suggest that the spines in the penumbra are basically the outward extension of the umbra. The spines and the penumbral filaments, 
together, are the basic elements that form a sunspot penumbra. 
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1. Introduction 


High-resolution observations of sunspots, which are the dark 
features on the solar surface, reveal that they contain small-scale. 


dynamically evolving structures, such as umbral dots ( 

Danielson 

1964 

[Riethmiiller et al. 2008[[2013[), light bridges dSo 

30tka et al. 

1997 

[Rimmele[20081 Shimizu[[2011[ [Lagg et al.[|20] 

4[), spines 

([Lites et al.|1993|), penumbral filaments (Tiwari et al. 

2013[ and 


references therein), and concentrated strong downflows at the 


spot’s periphery ( |van Noort et al.|2()T^ . In spite of the fast dy¬ 
namical evolution of these elementary features (~ 10 minutes 
to 4 h ours, see for example, |Sobotka|[l997t [Solanki & Ruedi| 
|2003| ), the sunspots are usually long-lasting, existing for a few 
days to several months, and only evolve slowly. To understand 
this relative stability, it is essential to study the global proper¬ 
ties of sunspots, in addition to exploring their fine-scale struc¬ 
ture. The global properties of sunspots are also of interest as 
constraints for global spot models, for studies of active region 
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magnetic fields and associated activity, for investigations of so¬ 
lar irradiance variations, and as proxies for starspots. 


In the past two decades, considerable progress has been 
made in understanding the three-dimensional (3D) sunspot 


structure in the solar atmosphere, both theoretically ([Rempel 

et al. 

2009a|bl [Rempel & Schlichenmaier[[2011[ Rempel[[2012) 

and observationally ([Solanki et al.| 19921 Lites et al.|1993 

Title 

et al. 

1993[[Stanchfield et al.[1997[[Westendorp Plaza et al. 

00 

a^ 

2001 

3|al[Solanki|2003l [Mathew et al.|2003l[Borrero & Ichimoto[ 

2011 

). Many of the depth-dependent global properties of the 


magnetic, thermal, and velocity fields in sunspots have been 
verified by those researchers. For example, within sunspots the 
magnetic field strength increases, whereas field inclination gen¬ 
erally decreases with depth, the temperature increases with both 
depth and radius, and the line-of-sight (LOS) velocity shows net 
upflows and downflows in the inner and outer penumbra, respec¬ 
tively. 


When looked at in detail, the more vertical fields in the 
penumbra (spines, according to [Lites et aL 1993] ) were found 
to be stronger than their surroundings. However, their thermal 
structure remains controversial. IWiehrl ( |2000l i associated dark 


1 



















































Tiwari, et al.: Global properties of a sunspot 


penumbral features with the stronger and more horizontal fields. 
Like |Westendorp Plaza et al.| ( |2Q01b| ), |Langhans et al.| ( |2005| ) 
found the str onger penumbral magnetic fi eld to be more vertical 
and brighter. Borrero & Ichimot^ ( 201l| ) found more horizon¬ 
tal fields (intra-spines) to be weaker and brighter in the inner 
penumbra and darker in the outer penumbra (see |Solanki|2003 


In spite of the large number of works on global sunspot prop- 


for a detailed review of earlier literat ure). These kinds o f con¬ 
troversies were recently addressed by [Tiwari et al. ( 2013| ), who 
show that the heads of sunspot penumbral filaments are brighter 
and contain vertical magnetic field, whereas their tails are darker 
and contain oppositely directed, strong, vertical magnetic fields. 
Only the darker regions that are less inclined to the vertical fields 
of the same polarity as the umbra, can be associated with the 
spines. It has also been confirmed that the magnetic field in the 
spines wraps around th e intra-spines in higher layer s, as antici¬ 
pated by the model s of jSolanki & Montavon (|1993|) and Spruit 


|& Scharrn^ ( |2006| ), and reported by jBorrero et al.| ( |2008| ). 


starts well wit hin the sunspot penumb ra and continues outside 
of it, whereas [Mathew et al. j ([2003 1) do not fin d a magnetic 


canopy structure anywhere. [Balthasar & Gomor^ ( [2008[ ) find the 


canopy only outside the penumbral boundary, in contrast to the 
above-mentioned researchers, but in agreement with, for exam¬ 
ple, [Giovanelfi] (p^80| ; [Giovanelir&Jon^ i 
( [1992[[1994[[1999[ ) and Wams et al.[ ( [199T 


1982); Solanki et al. 


Similarly, the relationship between the temperature and mag¬ 
netic field strength over the full sunspot is not fully under¬ 
stood ( [Solanki[[2Q03[ ). A non-linear relationship between tem- 
perature (and/or intensity) and magnetic fi eld is observed by 
Solanki et al. ( [1993[ ); [Mathew et al.[ ( [2004] ), which agrees with 


the anti-corr e lation found by other researcher s, e. g. |Kopp & 


Rabin[(p^9^ ; [Martmez Fillet & Vazque^ ( [1993[ ) and [Stanchfield 


et al.[ ( [T997f . The above authors attributed the non-linearity in 


Westendorp Plaza et al.[2001b|a 

[Solanki[2003MMathew et al. 

2003I Tritschler et al.||2004l |Bd 

lot Rubio et al.[|2004l [Borrero 


elementary, small-scale structures, such as umbral dots, spines, 
and penumbral filaments, has not been presented to date. 

Recent advances in the inversion of the high-quality 
spectro-polarimetric observations from Solar Optical 
Telescope/Spectropolarimeter (SOT/SP) provide us with 
an opportunity to take a fresh look at the above-mentioned 
issues in detail. In the present paper, we analyse the results of a 
depth-dependent inversion of a full sunspot and its surrounding 
region. T he same inverted re sults are analysed, as was done 
earlier by [Tiwari et al. ( 2013[ ). While we concentrated mainly 


However, many issues remain to be resolved. As an exam¬ 
ple, the magnetic canopy structure, which is a result of the ex¬ 
pansion of the magneti c field, has not been unambiguously ob- 
served and unde rstood. [Westendorp Plaza et al.[ ( [2Q01b[ ); [Rezaei[ 
et al. ( 2006[ ) and Borrero & Ichimoto[ ( [201 1[ ) find that the canopy 


on the internal structure of penumbral filaments in that paper, 
here we focus on the global and partly local properties of the 
inverted sunspot. 

This paper is organised as follows: Details of the observa¬ 
tions and the inversion process are given in the next section. In 
Section 3, we present the depth-dependent global thermal, ve¬ 
locity, and magnetic properties of the sunspot. The mutual de¬ 
pendence of physical parameters and inferred substructures of 
the sunspot are given in Section 4. In Section 5, we discuss our 
results, and finally, in Section 6, we present our conclusions. 


2. Observations and inversion 

We used a high spatial resolution spectro-polarimetric scan of 
the leading sunspot of the active region (AR), NOAA 10933, 
recorded with the SOT/SP ([Tsuneta et al.|2008||Suematsu et al.[ 


Hinode satellite (Kosugi et a 


[2008[ [Ichimoto et al.[|200^ Shimizu et al.[[2068[ ) onboard the 


.[2007 1. The positive polarity spot 


the (penumbral) lower temperature part of their scatter plots 
to the outer penumbral features. However, the different parts 
of the non-linearity in the scatter plots between field strength 
B and temperature T, e.g. a nearly constant temperature (of 
-5800-5900 K) for a range of magnetic field strength (0-2500 
G) found by [Westendorp Plaza et al.[ ( [2001b[ ), remains unex¬ 
plained. 

The magnetic field strength and inclination are anti¬ 
correlated, with more vertical fields having higher field strength 
(e.g. [Stanchfield et al.[[l997t [Westendorp Plaza et al.[[200Tb| ). 
This relationship mirrors the large-scale structure of sunspots. 
On smaller scales we expect a relationship that is affected by the 
local magnetoconvective structure, with stro ng opposite polarity 
fields at the tails of the penumbra l filaments ([Tiwa r i et al.[2013)); 
see also [Scharmer et al.[ ([2013[); jRuiz Cobo & Asensio Ramos[ 
( |20T3l ). 

The source of the Evershed flow ( [Evershed|1909[ ) is near the 
heads of the penumbral filaments, with the gas cooling down as it 
travels the bulk of the filaments and sink s as cool material at the 
tails of the filaments ( [Tiwari et al.|2013[ ). Based on this, a clear, 
systematic relationship between temperature and LOS velocity 
was found by [Tiwari et al.[ ( [20T^ for a standard penumbral fila¬ 
ment, where the upflows are hotter and the downflows are cooler. 
As an extension, in the present work, we perform this study on 
a full sunspot to investigate how this pattern, together with the 
cool spines, affects the structure of the full sunspot. 


of the NOAA AR 10933 that we use in the present study was 
observed close to the solar disk centre (/i=0.99) on 5 January 
2007 from 12:36 to 13:10 UT. The atomic lines contained in the 
data are the Ee I lines at 6301.5 and 6302.5 A. The SP scans are 
acquired in normal mode of SOT, with a spatial sampling of 0.16 
"pixel" ^ 

We perform the calibration of the Hinode (SOT/SP) data by 
using the standard SP_PREP routine, which is available in the 
SOLAR software package. The SP_PREP routine computes the 
thermal shifts in the spectral and slit dimensions first, and then 
applies the drift corrections for calibrating the data from Levels 
0 to 1 ( [Lites & Ichimoto|2013[ ). 

The observed continuum intensity map of the inverted re¬ 
gion is shown in Pig. The spot has an average diameter 
of about 50", and different aspects of it have already been 


studied in a number of publications: see for example, [Kubo 

[et al. ([2008); Tiwari et al.[ ([2009[); [Pranz & Schlichenmaier 

([2009); 

i Venkatakrishnan & Tiwari ([2009[[2010); 

Katsukawa & 

[Jurcak[ 

(201 0[);[Pranz[([20llj;[Borrero & Ichimoto 1 

2011) 

; Tiwari 

I20T21: 

i Riethmiiller et al. (2013[); van Noort et al. 1 

2013) 

; Tiwari 

[et al. (^ 

>013). 


2.1. Inversion 


We employ the output of the SPINOR inversion code ( [Erutiger 
et al.[[2000t [Erutig^[2000[) in the spatially coupled mode ([van 


Noort|201^[van Noort et al.|2013[ ) used by Tiwari et ST] ( [26l3[ ) 


to invert the calibrated Stokes profiles. The SPINOR code builds 
on the STOPRO line synthesis routines by [Solanki[ ( [1987[ ) and 
solves the radiative transfer equations numerically, under the as- 
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NOAA AR 10933 [05Jan2007:1236-131 OUT] 



0 10 20 30 40 50 60 

X [arcsec] 


Fig. 1. Continuum intensity map of the sunspot NOAA AR 
10933, observed by Hinode in normal scan mode of SOT/SP on 
5 January 2007 between 12:36 and 13:10 UT. The arrow indi¬ 
cates the direction of the solar disk centre. 


sumption of local thermodynamical equilibrium (LTE). The spa¬ 
tially coupled version of the SPINOR inversion code employs 
a modified Levenberg-Marquardt (L-M) algorithm for the opti¬ 
misation of a simplified, height-dependent atmospheric model, 
simultaneously over an extended field of view, taking the spa¬ 
tial degradation caused by the telescope into account. The re¬ 
quired spectro-polarimetric derivatives, with respect to the fitted 
atmospheric structure required by the L-M algorithm, are calcu- 
lated using response functions (j Landi DeglTnnocenti & Landij 
|DeglTnnocenti||1977} jRuiz Cobo & del Toro Iniesta||1992| ). To 


resolve the maximum amount of spatial structure, a Fourier in¬ 
terpolation of the data was made to allow the inversion to be 
carried out at an incre ases spatial sampling by a factor of two 
Ivan Noort et al.|2()T3] ). Thus, the pixel size for resultant physi¬ 


cal parameters used in the study is 0.08'' 

The vertical structure of the temperature T, magnetic field 
strength B, inclination y, azimuth LOS velocity vlos^ and 
microturbulent velocity Vmic are approximated using bicubic 
splines, controlled by three nodes, placed at log(T) = -2.5, -0.9, 
and 0, where r is the continuum optical depth at T = 6302.5 A. 

To help with the convergence of the solution, for every ten 
iterations, we re-initialised all slowly converging pixels by us¬ 
ing their nearest neighbours, and the whole solution was then 
smoothed over by convolution with a Gaussian smoothing func¬ 
tion. This significantly speeds up convergence if appropriate se¬ 
lection criteria for re-initialization are used. Although this pro¬ 
cedure can create a preference for a spatially smooth solution, 
it clearly results in a significant increase in the fit quality, i.e. a 
reduction in the difference between square of the data and the 
fitted profiles, i.e. a significantly lowerWhile, as in all inver¬ 
sions, the adopted solution may or may not represent the global 
optimum, it does appear to represent a workable optimum, since 
additional attempts to disturb and reconverge the solution did 
not significantly increase the accuracy of the fit, and resulted in 
a very similar solution. 

Examples of observed and fitted Stokes profiles for a typ¬ 
ical pixel, each one sample taken from penumbra and umbra. 


are shown in Fig[^ The atomic parameters used in the inversion 
are 6302.4936: log(gf)=-1.203, E/=3.6866 eV, and 6301.5012: 
log(gf)=-0.745, E/=3.6539 eV. Collisional broadening was cal¬ 
culated using the ABO tables ( [Barklem et al. 1998 ) and Fei 
abundance was set to 7.50. No molecular lines are included in the 
inversion since it makes the inversion process extremely slow. 
This sunspot was not particularly cool, meaning the inversion 
does not seem to have been severely affected by contamination 
through molecular lines (see, e.g. the fits for one of the darkest 
pixels in the sunspot umbra in Figj2]f _ 

A test was performed by [van Noort et al. ( [2013| ) to estab¬ 
lish the robustness of the inversion. This involved changing the 
node positions and adding different levels of noise in the ob¬ 
served and convolved profiles and resulted in very similar in¬ 
verted quantities, which suggests t h at the inversion result i s re¬ 
liable. Moreover, [van Noo^ ( |2012| ); [van Noort et al.| ( |2013| ) find 
good agreement between the inversion results and 3D radiation 
magnetohydrodynamics (MHD) simulations with the MURAM 
code ( |Vogler et al.|2005| ) for the velocities and fields, providing 
further support to the reliability of our inversion process. This 
result is supported by the fact that, for most pixels, the response 
functions are non-zero for all three chosen nodes, particularly 
for the one at log(T) = 0. 

However, establishing the error in the fitted atmospheric pa¬ 
rameters is a tedious task because of the spatial coupling of the 
parameters. A change introduced in a parameter at one pixel will 
cause the code to adjust the fit of that and other parameters in the 
neighbouring pixels. Therefore, determining the error in the fit¬ 
ted parameters would require a perturbation of each parameter 
at every pixel, followed by a full reconversion of the inversion 
in all the other pixels. This would demand a prohibitively large 
amount of computational power. 

The wavelength calibration of the data was carried out by 
requiring that the umbra of the sunspot, excluding the umbral 
dots, be on average at rest at r = 1. This resulted in a small 
systematic correction of the LOS velocity, vlos, of 100 ms“\ 
which was subtracted from the inversion results presented here. 

The 180° azimuthal ambiguity was res olved after the inver- 
sion by using the minimum energy method ( |Metcalf|1994t |Leka 


et al.|200^ . In Fig.|^ we display full maps of inversion results of 


some of the physical parameters (fitted temperature T, magnetic 
field strength B, field inclination y, and LOS vlos) at log(T) = 0. 
The physical parameters look smooth, in addition to having high 
contrast (e.g. compare the observed continuum intensity from 
Fig. with the one obtained from the inversion shown in Fig.[^ 
and discussed later). At log(T) = -0.9, where the response of 
the spectral lines 6301.5 and 6302.5 A to most of the physical 
parameters is largest, the physical parameters become smoother 
than at log(T) = 0. 


3. Global properties 

To investigate the global properties of the sunspot, we selected 
several contours of magnetic field strength of the spot, which 
was first smoothed by a factor of about twenty (20x20 pixels) 
to get smooth contours. Outer contours required smoothing by 
even larger factors (up to 50) because there are more fiuctuations 
along the azimuthal path that is present there. Ten contours out 
of a total of 21 are shown in red in Fig. The inner and outer 
penumbral boundaries are determined on the basis of intensity 
and are indicated by the yellow contours in Fig.|^ 

We note that both the outermost red contour and the outer 
yellow one differ quite significantly in some places. This is 
partly due to the different amount of smoothing applied to the 
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Fig. 2. Examples of observed (black) and best-fit (red) Stokes profiles: I, Q, U, and V. Upper row displays Stokes profiles for a 
typical pixel in the sunspot penumbra (location: X48", Y41.6"). Lower row contains Stokes profiles for one of the darkest pixels in 
the sunspot umbra (location: X36",Y32"). 
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Fig. 3. Maps of selected physical parameters (at log(T) = 0) obtained from the inversion: (a) temperature T, (b) magnetic field 
strength B, (c) magnetic field inclination y, (d) line-of-sight velocity vlos- The arrow in the temperature map points to the solar disk 
centre. Except for field inclination y, thresholds have been used for minimum and maximum values of the colour bars. The extreme 
values returned by the inversion are 3853 to 7173 K for T, 10 to 4000 G for B, and -11.5 to 18.7 kms"^ for vlos- 
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Fig. 4. Ten of the 21 contours of magnetic field strength over¬ 
plotted on the inverted continuum intensity map. Red contours 
are based on the smoothed magnetic field strength. Yellow con¬ 
tours are based on the intensity and represent the boundaries of 
the sunspot penumbra. The greyscale bar ranges from 0.1 to 1.3; 
the extreme values of IdIc<Quiet Sun> returned from the inver¬ 
sion are 0.097 and 1.468. The black arrow on the upper right part 
of the image points towards the solar disk centre. 


field strength and continuum intensity prior to determining the 
contours, but partly also to a mismatch between field strength 
and intensity near the penumbral boundary. We looked care¬ 
fully at these regions and found that many of them are the tails 
of penumbral filaments with high-speed downflows, and con¬ 
tain very strong fields with a polarity opposite to that of the 
sunspot umbra ( jvan Noort et al.|20T^ . In most of these regions 
the tails of several penumbral filaments seem to converge. The 
most prominent excursion of the outermost red contour beyond 
the yellow contour (X 10", Y 45") is in the direction of a 
group of pores of the same polarity, i.e. in a direction in which 
the penumbra is strongly distorted. 

3.1. Thermal properties 

To explore the variation in the continuum intensity and temper¬ 
ature along the radius of the sunspot, we averaged these param¬ 
eters azimuthally along the given contours (i.e. at the three opti¬ 
cal depth positions for temperature that correspond to the nodes 
used for the vertical spline interpolation of the atmosphere). The 
mean values, along with the root mean square (RMS) variation 
of the sample along the path of the contours (not to be confused 
with the standard error to the mean), are plotted in Fig. The 
RMS variations are quite significant, indicating the presence of 
fine-scale structures along the path of the contours. The contours 
are based on the magnetic field strength and are not equidis¬ 
tant. Because of this, the corresponding error bars are also not 
equidistant. 

The brightness of the umbra increases with radial distance 
from the centre of the spot. We find that the umbra has an az¬ 
imuthally averaged intensity of 20-40% of the quiet Sun in the 
continuum at 6302.0 A, in agreement with earlier findings (e.g. 



Normalised Radius 


Fig. 5. Upper panel: Radial distribution of the azimuthally av¬ 
eraged continuum intensity normalised to averaged quiet Sun. 
The error bars represent the RMS variation along the azimuthal 
path of the contours. Lower panel: Radial dependence of the az¬ 
imuthal averages of temperature at three optical depth positions. 
The two vertical dashed lines represent average radial positions 
of the boundaries of the sunspot penumbra. 


|Westendorp Plaza et al.|[^01b| ). The darkest part of the um¬ 
bra has continuum intensities as low as 10% of the quiet Sun. 
The penumbra has a brightness of 45-85% of that of the av¬ 
eraged quiet Sun. However, these values depend on the crite¬ 
rion used to determine the umbra -penumbra boundary and on 
the wavelength. (See, for example, [Mathew et al.| ( |2Q03 ) for in¬ 
frared wavelengths). They also depend on the size of sunspots 
( [Mathew et al.| 2007 1 and references therein), whereas the pur- 
ported dependence on the phase of solar cycle remains inconclu- 


sive (jAlbregtsen & Maltby||1978 

|1981[ Albregtsen et al.||1984l 

Maltby et al.|1986[ Norton & Gil 

man|2004[ [Penn & Livingston 

20061 IMathew et ai|2007|l. 

hally averaged temperature at 
vn in the lower panel of Fig. ^ 
uum intensity, displayed in the 
This is to be expected for an 

The behaviour of the azimut 
the deepest node log(T) = 0, sho\ 
closely mimics that of the contin 
upper panel of the same figure. 


LTE inversion result. The temperature drops more rapidly from 
penumbra to umbra in the deeper layers than higher up, suggest¬ 
ing a lower vertical temperature gradient in the spot. Whereas in 
the penumbra this vertical temperature gradient is only slightly 
weaker than that in the quiet Sun, it is much weaker in the umbra. 
The latter result is in good agreement with the literature (see, 
e.g. |Schroter||197ll|Maltby et al.|[T^ l^ados et al.|[T^ 


5 














































































































Tiwari, et al.: Global properties of a sunspot 


Westendorp Plaza et al.||200fbt |Solanki|[2003[ [Tritschler et al. 


2004| ). 


The temperature shows the largest fluctuations at log(T) = 
0. Since the temperature is very well constrained in this layer, 
largely by the continuum intensity, such strong fluctuations are 
an intrinsic property of the sunspot atmosphere, and are caused 
by the flne-scale structures in it. 


3.2. Velocity stratification 


Our sunspot was located almost at the centre of the solar disk 
(within 5°), allowing us to treat the LOS velocity as the verti¬ 
cal velocity, and to ignore projection effects. We And that the 
deepest layer in the inner penumbra is dominated by upflows 
(see Fig. [^), in agreement with earlier findings. This is as ex¬ 
pected from the dominating presence of heads of penumbral fil¬ 
aments, which are found in this location of the spot, and which 
contain strong upflows ( [Tiwari et al.|2013 ). At the outer bound¬ 
ary of the sunspot penumbra, strong downflows dominate in 
the deepest layer, again as expected from the known concen¬ 
trations of the tails of penumbral filaments with strong down¬ 
flows that are found at this location ( [van Noort et al.j[20131 


Tiwari et al.|2013[). This large-scal e radial structure of flows 


was 


seen earlier by 

Rimmele ([1995a[); jSchlicher 

imaier & SchmicW 

chler et al.|p004|); 

o 

o 

o 

(N 

Westenc 

iorp Plaza et al.[p001a|)nTrits 

Rimme 

e & Marino[ 

([2006J; Ichimoto et al. 

P007a); Sanchez[ 



Normalised Radius 

Fig. 6. Azimuthal averages of vlos at three optical depths plot¬ 
ted vs normalised radial distance from centre of the spot (solid 
lines). The RMS variations along each azimuthal contour are in¬ 
dicated for log(T) = 0 and -2.5 by the red and blue vertical bars, 
respectively. The two vertical dashed lines have the same mean¬ 
ing as in Fig.[^ 


In Fig. we present the azimuthally averaged vlos at the 
three depth node positions vs radius. This plot clearly shows the 
domination of upflows in the inner penumbra with the average 
velocity in the deepest layer around -300 ms“^ and the even 
greater domination of downflows in the outer penumbra, where 
the average velocity reaches values as high as 1300 ms“^ just 
outside the sunspot boundary. The rapid increase in the RMS 
variations with depth is in good agreement with the results of 
other studies of the flne structure of penumbrae (see for example 


[Tiwari et al.|2013[ and references therein). The increase in RMS 
fluctuations is nearly linear, with radial distance up to the outer 
boundary of the spot, which indicates increased inhomogeneity 
in the vertical velocity. 

The average downflow in the deepest layer increases in 
strength up to a normalised radius of ~ 1.05, outside the vis¬ 
ible boundary of the sunspot. The increase of the velocity be¬ 
yond the boundary of the spot is a result of some of the strongest 
downflowing regions not being inside the contour that marks the 
sunspot boundary, as pointed out in the beginning of Section 
Furthermore, the downflows continue well beyond the boundary 
of the sunspot, in agreement with Boemer & Kneer[ ( [1992[ ). 

The average velocity in the umbra is zero at log(T) = 0, w hich 
was set by the velocity calibration described in Section [XT] The 
average velocity in the higher layers of the umbra remains close 
to (but different from) zero throughout the sunspot. In partic¬ 
ular, there is a systematic weak downflow in the umbra above 
log(T) = 0 , wh ich increases with height, although, as we discuss 
in Section [5T| the reliability of this result remains uncertain. 

The average velocity in the penumbra is a significant down¬ 
flow of about 350 m s“^ at log(T) = 0, as expected from the pres¬ 
ence of strong downflows at the outer boundary of the penumbra. 
This average downflow decreases with height, reaching zero at 
log(T) = -2.5. 


3.3. Magnetic properties 

3.3.1. Azimuthal averages of field strength B and inclination 

r 

Azimuthally averaged values of B and y, along with their RMS 
variations along the contours, are plotted in Fig. [T] The large 
RMS variations are again indicative of the presence of flne-scale 
structures along the path of the contours, although in the um¬ 
bra the RMS variations tend to indicate the larger-scale inhomo¬ 
geneity of B and y. The lower values of the error bars go below 
zero for the held strength for the contours outside the penumbral 
boundary because the distribution of the held strength along the 
path of these contours departs significantly from Gaussian. 

We can see that the decrease in the average magnetic held 
strength is nearly linear with radial distance from the centre of 
the spot. The drop in held strength at log(T) = 0, from 2800 G (in 
umbra) to 700 G (at outer penumbral boundary), is slightly larger 
than at the other two heights. Also noticeable is the substantial 
weakening of the held strength in the two lowest nodes outside 
the sunspot boundary. 

The general trend toward a decrease in the held strength with 


Plaza et al. pOOlbJ; [Mathew et al.[ P003[), and [Borrero & 

Ichimoto 

(201 1[), but the height dependence of B in our plots, as 

discussec 

later, does not match those of 

Westendorp Plaza et al.[ 

(|2001b|lJ 

Borrero & Ichimoto [POIIJ, or 

Mathew et al.[p003[). 


averaged magnetic held strength increases with height: 1) just 
outside the visible boundary of the sunspot penumbra, and 2) in 
the region between radial positions 0.43 and 0.63. These are the 
locations of canopy structure and positive he ld grad ient, respec¬ 
tively, which we describe in the next section ( 3.3.2[ ). 


The azimuthally averaged held inclination increases with 
radius, in agreement with the behaviour of an expanding flux 
tube. The averaged magnetic held inclination consistently ex¬ 
hibits more vertical fields in the higher layers in the umbra and 
penumbra. The magnitude of this vertical gradient, e.g. 5-10° 
between the lower and middle photosphere, is sizable. Outside 
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Fig. 7. Azimuthal averages of B (upper panel) and y (lower 
panel) along the sunspot radius at three optical depths. Also 
shown are RMS variations for the uppermost and lowermost 
node positions (blue and red error bars, respectively). The two 
vertical dashed lines represent the inner and outer boundaries of 
the sunspot penumbra (on average). The two dotted lines outline 
the radius between which a positive field gradient is seen (i.e., 
field strength increasing upward). 


the sunspot, however, the field in the deepest layer becomes 
more vertical and azimuthal averages of inclinations at all three 
heights are nearly the same. 

The fiuctuations in both the field strength and inclination are 
largest in the deepest layer, and those in the inclination increase 
towards the outer boundary of the sunspot. 


3.3.2. Field gradients and canopy structure 

To estimate the magnetic field gradient, we need to obtain the 
field strength on a geometrical height (z-) scale, whereas the in¬ 
version code delivers the stratification of each physical parame¬ 
ter on a continuum optical depth (r-) scale. The transformation 
of optical depth scale to geometrical height is not straightfor¬ 
ward because of the unknown optical corrugation of the surface 
over sunspots. We therefore carried out this transformation by 
assuming hydrostatic equilibrium. We are well aware that this is 
a gross simplificati on since it neglects the in fluence of the mag¬ 
netic field (see also |Puschmann et ar| ( |2010| )). 

The average vertical field gradient dBjdz in the sunspot um¬ 
bra near the deepest node, where the highest field strength is 
about 4 000 G, was determined to be -1.4 G km“^. The gradi¬ 
ent becomes weaker with height and reaches a value of -0.95 G 
km“^ at log(T) = -2.5. The magnitude of the field gradient found 



0 10 20 30 40 50 60 

X [arcsec] 


Fig. 8. Gradient of magnetic field computed between the geo¬ 
metrical height corresponding to the two deepest nodes (log(T) 
= -0.9, 0). A positive field gradient represents field strength in¬ 
creasing with height. 


for our sunspot is equivalent to that found in 3D magnetohydro¬ 
dynamic (MHD) simulations (-1.5 G km“^ for a sunspot whose 
strongest magnetic field is close to 4 000 G at r = 1: M. Rempel, 
2013, private communication). 

The azimuthally averaged magnetic field shows a positive 
gradient (opposite to that seen in the umbra) in the inner-middle 
part of the penumbra (roughly between radius = 0.43 and 0.63), 
and just outside the sunspot (see upper panel of Fig. [7]). The pos¬ 
itive gradient observed outside the sunspot penumbra can clearly 


be associated with the magnetic canopy structure ( 

Giovanelli 

T^IGiovanelli & Jones|1982l|Solan]d et al.|1992l 

Solanki & 

Schmidt||1993[ Solanki et al.||1994[ 1 19991), whereas the inverse 


gradient seen in the inner-middle part of the penumbra’s deep¬ 
est layers must have a different source. Figure clearly reveals 
the locations of the positive field strength gradient present in the 
inner-middle penumbra and at the outer penumbral boundary. 
Such an inverse gradient is a lso seen in Vacuum Tower Telescope 
(VTT: [Schroeter et al.|1985| ) data by |Joshi et al.| ( |2015| , and Joshi 


(2014), who analysed this particular property of the sunspots in 
detail. They also found similar positive gradients in the inner 
penumbrae of MHD simulations of sunspots. 

Although, locally, a positive field gradient is clearly visible 
between the optical depth layers log(T) = -0.9, and 0, this is not 
the case for the azimuthal averages outside the spot, where the 
negative gradient contributions dominate the positive ones. Only 
in the upper node is the average field strength greater than in the 
deepest layer (see upper panel of Fig.[^. 

The strongest negative field gradients concentrated locally 
near the outer spot boundary (see Fig. are associated with 
strong downflows, having the opposite polarity to the umbra. 
The very strong field strength in these regions in the deepest 
magnetic layers (some of them reaching up to 4 kG, the value 
set a priori as an upper limit in the inversion) can probably ex¬ 
plain the strong n egative field gradients seen in these locations 
(see, for example, van Noort et al.|20f3 for details). 
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3.3.3. Field azimuth and twist in the sunspot 

A field azimuth map, obtained at the middle node (log r = -0.9), 
is shown in the left-hand panel of Fig. It displays a generally, 
radially directed field with a considerable amount of small-scale 
structure. 


% [°] 
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Fig. 9. Left panel: field azimuth map at log(T) = -0.9. Right 
panel: twist map calculated from the difference between the 
measured field azimuth and radial directions from the sunspot 
centre. To avoid contributions of the azimuth from noisy signals 
in the quiet surroundings, only pixels lying within the sunspot 
penumbral boundary have been used to compute the twist. 


To estimate the magnetic twist (average deviation of the field 
azimuth from the radial direction), we have used only pixels ly¬ 
ing within the sunspot boundary. Twist is plotted in the right- 
hand panel of Fig.|^ We excluded the quiet-Sun region to avoid 
the influence of noise and fluctuations. The average twist value 
of the full sunspot in the b ottom layer comes out at -3.8°. This 
value for twist agrees with [Tiwari et al.| ( [2009| ), who estimated a 
twist of the same spot a day before by using a potential field as 
reference, instead of the radial field. Also, in agreement with 
these findings, it is clear from the right-hand panel of Fig. 
that much larger and oppositely directed twists dominate locally. 
Sometimes the local twist reaches a value of |A;^| > 90°. The es¬ 
timated azimuthally averaged twist with radius at the three node 
positions are plotted in Fig. However, the large twist in the 
umbra seen in Fig.j^is partly due to the geometric and magnetic 
centres of the spot not being completely co-located. For this rea¬ 
son only, the twist of the penumbral field is plotted in Fig.p^ 

In the penumbra, there is a trend toward increasing twist with 
both radius and height. In the middle penumbra, the average un¬ 
signed twist in the lower two nodes is generally less than 5°. 
Although the variations in the field azimuth along the contours 
are large, the estimated error in the mean value is not, owing to 
the large number of pixels along the contours. The error in the 
mean is indicated in Fig.[T^and is clearly much smaller than the 
recovered twist values. However, the error bars are not important 
for the contours if they have lots of pixels with |A;^| > 90°. 

It is also worth mentioning that the field azimuth is the most 
fluctuating parameter obtained from the inversion. It is also the 
one most strongly affected by noise. Nonetheless, Fig. 9 does 
indicate that there could be a gradient in the twist of the sunspot 
field, both with radial direction and with height. 



Fig. 10. Radial distribution of twist, the azimuthally averaged 
difference A;^, for the sunspot penumbra at three nodes. All ver¬ 
tical lines are the same as in Fig.[^ The RMS variations are very 
large along the azimuthal path and not shown here. The error 
bars represent the standard error of the mean along an azimuthal 
contour. 


3.4. Mass flux 


For a reliable calculation of the mass flux, we require that the 
density and the LOS velocity be on a geometrical height scale. 
Not only can the relative geometric height differences between 
different points in the field of view not be inferred by the inver¬ 
sion code, the density and vertical height scale is also calculated 
by assuming hydrostatic equilibrium, an assumption that is often 
not valid in a highly magnetised structure. 

When bearing the above limitations in mind, and assuming 
that T unity lies at the same z at all positions in penumbra, a 
crude estimate of the mass flux over the full sunspot gives about 
2.5 times more downfiowing mass than up flowing, which agrees 
with the estimates of jWestendorp Plaza et al. (|1997| [200laj ) for 
a sunspot and with the results of [Tiwari et al.| ( 2Q13| ) for an aver¬ 
aged sunspot penumbral filament. The amount of excess down¬ 
flows decreases with height, but even at log(T) = -2.5, the ex¬ 
cess is a factor of about 1.3. 


Although a partial contribution to the downfiowing mass 
from the inverse Evershed flow effect through spines cannot be 
ruled out, this result suggests that the assumptions made when 
deducing the mass flux, e.g. hydrostatic equilibrium, log(T) = 0 
located at the same z for all penumbral pixels, cannot be true. 
Therefore, a true geometrical height computation is necessary 
for estimating the mass flux over a sunspot reliably. 

Finally, the net mass flux depends on the zero level of the 
velocity. We set the velocity at the lowest node in the dark parts 
of the umbra to zero, since the density is highest at that layer. 
However, if we force the velocity to be zero at the central node 
(log T = -0.9), then the mass flux excess at log(T) = 0 reduces 
to 1.7 (down from a factor of 2.5). 


4. Mutual dependence of physical parameters and 
sunspot substructures 

In this section, we plot two dimensional (2D) histograms of one 
physical parameter relative to another and investigate their mu¬ 
tual dependences. Then, by isolating different populations in the 
histograms, we identify specific substructures of the sunspot. 
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Fig. 11. Left panels: 2D scatter plots of 5 vs y for all points within the sunspot’s boundary at the three optical depth nodes (from 
top to bottom, log(T) = -2.5, -0.9,0). Right panels: the same but only for pixels in the penumbra. The vertical dashed line in each 
panel represents the location of y = 90°. A black box is drawn at log(T) = 0 to identify/outline the umbral pixels, spine pixels in (f) 
(see Fig.[T^a) and (b) for where these pixels are located in the sunspot). 


4.1. Field strength B versus inclination y 


The 2D histograms of B versus y of the full sunspot at the 
three height node positions are shown in the left-hand panels 
of Fig. [m A general anticorrelation between B an d y is found 
at all depths for y < 90°, in general agreement w ith|Stanchfield 




1997| ); |Westendorp Plaza et al.| ( |2001b| ) and |Mathew et al 


( |2004| ). A^eak positive correlation, not previously reported, is 
noticeable at all heights in the range y = 90° - 180°, in spite of 
the relative scarcity of points. 

The two peaks in the histograms at all the three heights may 
be interpreted as belonging to the sunspot umbra (the popula¬ 
tion in the black box shown for log(T) = 0 in Fig. [TT^) and the 
penumbra (the population near y = 90°). In an attempt to clarify 
this picture, we recreated the histograms but only for penum- 
bral pixels, shown in the right-hand panels of Fig.[TT] From the 
residual population of the black box, it is clear that the popu¬ 
lation, which we were relating to umbra, consists not only of 
umbra, but also, to a significant extent, of penumbral pixels. By 
marking all pixels contained by the black box in Fig. [^e), in 
red in Fig. 12 a), and in blue in Fig.p^b), it becomes clear that 


some of the pixels in the black box belong to the spines of the 
penumbra, which appear to be true extensions of the umbra into 
the penumbra. The two populations of umbral/spine pixels and 
penumbral filament pixels, which are quite distinct at log(T) = 0, 
increasingly merge with height, although the umbral/spine pix¬ 
els exhibit much less scatter at log(T) = -2.5 than the filament 
pixels. 


The remaining penumbral pixels, including the second pop¬ 
ulation characterised by nearly horizontal fields with an av¬ 
erage field strength of 1 kG, consist of penumbral filaments. 
Penumbral filaments, whose bulk contains a horizontal field with 
a strength of ~1 kG ( [Tiwari et al.|2013] ), clearly form the domi¬ 
nant part of the penumbra. The right halves of the frames, with 
y > 90°, are populated by the filament tails where the field is 
stronger, with a polarity opposite to that of the umbra. The op¬ 
posite polarity is most prominent at log(T) = 0 but is also visible 
in the higher nodes (albeit for fewer pixels), consistent with the 
field bending back down into the photosphere. 
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Fig. 12. (a) Map of magnetic field strength at log(T) = 0. The pixels with B > 1600 G & y < 50° are highlighted in red (i.e. the 
pixels in the box outlined in black in Figs.[^e) and 10(f)). (b) Field inclination map at log(T) = 0. with pixels satisfying the same 
conditions as applied to the field strength map in panel (a) highlighted in blue. The colour scheme for the field inclination map is 
chosen for consistency with the field inclination map depicted in Fig. (c) Same map as displayed in (b) but pixels with vlos > 
3.5 kms“^ are highlighted in black, revealing tails of penumbral filaments, (d) Map of temperature at log(T) = 0. Highlighted in 
white are pixels satisfying T > 6400 K & 1000 < 5 < 2500 G, basically representing heads and warm parts of bulks of penumbral 
filaments. The dash-dotted contour in each panel represents the umbra-penumbra boundary. 


4.2. Temperature T versus magnetic field strength B 


Biermann] ( |1941| ) was the first to realise that the presence 
of a strong mag netic field in a sunspot could be responsible 
for its darkness. Alfven ( |1943| ) predicted a relationship be¬ 
tween the magnetic field strength and temperature, which has 
been extensively investigated both theoretically by, for example, 
Outml ( p^63] ); [Plckil ( p^TO] ); |Cowling| ( [T9761 ); |Maltby| ( p^TT] ) 


and [Spruit et al.| ([1990|), and observationally b y, for example 
Chou|^1987|); Ma ^rtin ez Fillet & Vazque^(|1990|); Kopp & Rabin 


; [Martinez Fill et & Vazqu ez] (|1993| ); Westendorp Plaza 
|et al.|j2001b| ) and |]^thew et al.[^2004| ). However, the relation¬ 
ship between B and T of sunspots, particularly the correspon¬ 
dence between different parts of scatter plots found by these au- 
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thors and feat ures of sunspots, has not yet been understood (see 
Solanki|2Q03[ for a detailed review of the subject). 


The 2D histograms of T versus B in Fig. [fallow us to re¬ 
visit the thermal-magnetic relationship of sunspots. The temper¬ 
ature and its spread over different parts of the sunspot decrease 
rapidly with height. At first glance, the plots for the full sunspot 


in the left-hand panels ap 

3ear nearly identical to the results ob- 

tained by[Kopp & Rabin|( 

1992[);[Solanki et al. ([1993J; [Martinez 

Fillet & Vazquez[([1993[); 

Stanchfield et al. ([1997|);[Westendorp 

Plaza et al. [ ([200lb[);[Penn et al.[([2003[) and[Mathew et al.[([2004[). 


However, a closer look reveals noticeable differences, mainly in 
the right-hand parts of the plots, which are expected to be pop¬ 
ulated by penumbral pixels. To confirm this identification, we 
plotted the same histograms at the three heights but only for 
penumbral pixels, shown in the right-hand panels of Fig.The 
most densely populated area is at around 1 kG field strength and 
a temperature of 5 900 - 6 400 K at log(T) = 0. 


Although the general impression of Fig. 13 is of an anti¬ 
correlation (also in the penumbra, where the correlation coef¬ 
ficient is around -0.4), we do find a localised weak positive 
correlation between the two parameters, particularly visible in 
the region outlined in Fig. [^f). The positive correlation be¬ 
tween the temperature and the magnetic field strength is prob¬ 
ably caused by the brighter heads of the penumbral filaments, 
which are brighter than their surroundings and often contain a 
stronger magnetic field ( [Tiwari et al.|2013) . To confirm this, we 
identified the locations of the pixels lying in the black box, and 
highlighted them in Fig.[^d), which confirms that they are, to 
a significant extent, the locations of the heads of penumbral fila¬ 
ments; some pixels, e.g. from the bulk of filaments, also appear 
to be present. An additional scatter plot (not shown here) of B vs 
r, showing only highlighted pixels, confirms this weak positive 
correlation. However, not all the heads of penumbral filaments 
are highlighted by the selection criteria described in Fig.[T^d). 
For example, some heads may contain a temperature of 6400 K 
or less, which then fall outside the box towards lower tempera¬ 
tures and are not highlighted in Fig.fT^d). 

Another anomalous, rather scattered, population in the 
Fig.[^f) is formed by the pixels above 2500 G and 5500 K (in¬ 
dicated by an arrow). This is formed by the tails of penumbral 
filaments, which have stronger fields a nd are somewhat da rker 
than the heads of filaments. As shown by Tiwari et'n^ ( |2013| , the 
tails of filaments show an enhancement in temperature, some- 
times becoming h otter than parts of their bulk (see Fig. 5 of 
Tiwari et 3T]|2()13| ). The panels (e) and (f) of Fig. [T3] show that 
the strongest fields in the sunspot belong either to the umbra or 
to the tails of penumbral filaments. 

The highlighted pixels for tails of penumbral filaments in 


Fig. 12 c), by the condition vlos >3.5 kms \ as described 


later, also accommodate the pixels with T{> 5500^) and B{> 
2500G) of the population indicated by an arrow in Fig.[T3]^f). 


4.3. LOS velocity vi^os versus temperature T 

The upfiowing gas is expected to be hotter than the downfiow- 
ing gas if these fiows are driven by convection. Two dimen¬ 
sional histograms of vlos versus T for the full sunspot at the 
three node positions are displayed in the left-hand panels of 
Fig. The majority of the coolest pixels (mostly belonging 


to the umbra) are concentrated around zero velocity at all three 
heights (the horizontal branch of the histograms to the left), with 
some deviations in the deepest node because of umbral dots (e.g. 
Riethmuller et al.|2013| ) and faint light bridges (e.g. [Lagg et ak] 
2014| ) containing up- and downflows. 


The broad distribution of points above roughly 5500 K in 
all panels is the penumbral contribution, as confirmed by similar 
histograms of only the penumbral pixels and as depicted in the 
right-hand panels of Fig.[^ It is clearly visible that the upfiows 
at log(T) = 0 are, on average, hotter than the downfiows. The 
average temperatures of up- and downfiows are ~ 6300 and 6050 
K, respectively, for full penumbral pixels, and 6600 and 6200 K 
for the pixels outside the two solid horizontal lines. There is a 
weak tendency for the strongest downfiows to be warmer than 
the weak downfiows. 

The hot upfiows and cool downfiows, particularly the ex¬ 
tended populations below and above the horizontal lines, respec¬ 
tively, in the panels (e) and (f) of Fig. can readily be associ¬ 
ated with the heads and tails of penumbral filaments. To confirm 
that the extended population above the upper horizontal line in 
Fig. [^f) belongs to the tails of penumbral filaments, we high¬ 
lighted those pixels in Fig.[^c). Similarly, we tested if the pop¬ 
ulation below the lower horizontal line in Fig. [T^f) belongs to 
heads of penumbral filaments, and found that they are accommo¬ 
dated in the highlighted pixels in Fig.[^d). The pixels below the 
lower solid line in Fig.fT^f) are only a small subset of the pixels 
highlighted in Fig.[T^d). 

In Figure [^f), the spines are cooler and contain very low 
velocities. These spines get cooler with height and other penum¬ 
bral pixels merge with them. In the topmost layer, in Figure 
[Hb), unexplained weak upfiows can be seen that appear to be 
cooler than the downfiows. Such a phenomenon is well-known in 
the quiet Sun, and known as reversed granulation, with granules 
becoming cooler and intergranular lanes becoming hotter above 
log optical depth of roughly -2. We speculate that this (weak) 
reversal, seen at log(T) = -2.5 in the penumbra, is the signature 
of the same effect (described by, for example, [Nordlund et al. 
( |2009| ) and is the result of the adiabatic cooling of the expanding 
upfiowing gas). 


4.4. LOS velocity vi^os versus magnetic field strength B 

Figure shows the 2D histograms of vlos vs B for the whole 
spot (left-hand panels), and for the penumbra only (right-hand 
panels). The most prominent features in the histograms are the 
two branches towards the upper and lower right corners of the 
domain, most clearly visible at log(T) = 0. These are readily as¬ 
sociated with the tails and heads of penumbral filaments, respec¬ 
tively, and suggest that the strongest magnetic fields are not only 
found in the umbra, but also in the deepest layers of the fastest 
downfiows, which are expected to lie in the tails of penumbral 
filaments, based on the work of [van Noort et al.j ( |2013j ) and 
[Tiwari et~ST] ( |2013j ). 

To confirm the above prediction of the positions of pixels be¬ 
longing to the heads and tails of penumbral filaments, we over¬ 
plotted separately all the pixels outside the two solid horizon¬ 
tal lines shown in Fig.[n];f). As shown in Fig.[^c), the pixels 
with Vlos >3.5 kms^^create a good match with the tails of 
filaments. The pixels with vlos < “2.5 kms“^ are completely 
accommodated within the highlighted pixels of Fig.fT^d). 

There is a significant difference of 300 G (1 500 G vs 1 800 
G) in the field strengths of average upfiows (-3.9 kms“^) and 
downfiows (5.8 kms“^) computed outside the solid lines, the 
latter being stronger. This difference can generate an outf low via 
a siphon flow, if it refers to the same geometrical height ( [Meyer 
|& Schmidt|[1968[ [Montesinos & Thomas[[1997] ). However, the 
higher average temperature (6600 K) of upfiows than downfiows 
(i.e. 6200 K) supports the convection mechanism as the driver of 
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Fig. 13. Left panels: 2D scatter plots of .S vs L for all pixels. Right panels: The same but only for penumbral pixels inside the 
sunspot, shown at the optical depths of the three nodes. The black box in (f) outlines the heads of some of penumbral filaments (see, 
Fig.[^d) to find the spatial locations of these pixels). The arrow in (f) points to the population that belongs to tails of penumbral 
filaments that are subsets of highlighted pixels in Fig.p^c)). 


the Evershed flow. Thus, a conclusive statement about the driv¬ 
ing mechanism of the Evershed flow cannot be made. 

The reduced population of spine pixels in the left- and right- 
hand panels of Fig.[T^at all heights again shows that the sunspot 
umbra and spines exhibit similar physical properties. A rapid re¬ 
duction in both the up- and downfiows with increasing height 
is in agreement with the literature (see for example, [Westendo^ 
[Plaza et al.|20QlatrSolanki|2()03| and references therein). A weak 
tendency of a shift from overall dominant downfiows to weak 
dominant upfiows with increasing height can be observed, how¬ 
ever. 

5. Discussion 

We have performed depth-dependent, spatially coupled inver¬ 
sions of a disk-centered sunspot, observed by Hinode SOT/SP, 
and have presented the depth-stratified thermal, velocity, and 
magnetic atmospheric structure of it. We have also looked at 
how the small-scale structures fit in with the global behaviour 
of the sunspot and have introduced simple thresholds in individ¬ 
ual parameters to isolate some of the fundamental constituents 


of the sunspot penumbrae. In the following, we interpret our re¬ 
sults and discuss them in the context of earlier work available in 
the literature. 

5.1. Global properties 

The continuum intensity and temperature at logir) = 0 display a 
very similar trend, an obvious increase with radius starting from 
the centre of the spot. The horizontal gradient gets smoother for 
higher layers. Results imply that the vertical temperature gradi¬ 
ent in the sunspot penumbra is slightly smaller than that in the 
quiet Sun, but is much larger than that in the umbra. These re¬ 
sults are in good agreement with findings in the literature (see. 


for example, Schroter 

197 1[ 

Maltby et al.|1986[ [Collados et al. 

1987 

, Lites et al. 1993 

Stanc 

ifield et al.[1997[[Westendorp Plaza 

et al. 

2001bl|Solanki|2003l|Tritschler et al.|2004D. 


The azimuthally averaged magnetic field strength decreases 
with radial distance from the centre of the sunspot (2800 G) to 
the outer penumbral boundary (700 G at r unity level), in qual- 
itative agreement with earlier find i ngs, e.g. jWestendorp Plazaj 
jet al.| ( |2^Tb| ); |Mathew et al.| ( |2003| ); [Borrero & Ichimoto| ( |201 Ij ). 
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Fig. 14. Left panels: 2D scatter plots of T vs vlos for all points within the sunspot boundary at the three optical depth nodes. The 
horizontal dashed line in each panel indicates the zero velocity level. In each panel, the points at lower temperature belong to the 
sunspot umbra and also to spines, while those at higher temperature belong only to the penumbra. Right panels: the same, but only 
showing penumbral pixels. The two horizontal solid lines are used to separate, in an approximative manner, the heads (lying below 
the lower solid line) and tails of penumbral filaments (above the upper solid line) from the other points. The dashed vertical line in 
(f) represents T = 6400 K, a condition applied to isolate the heads of penumbral filaments in Fig. [^d), although it also includes 
points from the bulks of the filaments (points to the right of the vertical dashed line lying between the two horizontal solid lines). 


However, some quantitative differences can be seen, beyond the 
well-k nown difference in the maximum B of spots (e.g. |Schad| 
|2013| ), which is known to depend on the size of spots. For exam¬ 
ple, the averaged field strength values of about 2300 G and 2500 
G for umbra, an d 500 G and 700 G for outer penu mbr al bound¬ 
aries, found by [Westendorp Plaza et al. ( |2001b| ), and Mathew] 
et al.lpOO^] ), respectively, might be due to the variation of prop¬ 


et al. 


erties from one sunspot to another (see also |Solanki ( |2003| )) but 
may also reflect differences in the spectral lines used and the 
spatial resolution of the data. The maximum B in our sunspot is 
comparable to that of the sunspots with similar areas in the work 
of |Sctmd| ( |MT3] ). 

The vertical gradient in the photospheric layers of the umbra 
of -1.4 G km“^ is comparable to that obtained in the most recent 
MHD simulations (e.g. -1.5 G km“^: M. Rempel, 2013, private 
communication). In the literature, vertical fi eld gradient values 
in sunspot umbrae varying from -4 G km“^ (Westendorp Plaza 


to -1.5 G km ^ (Collados et al.|1994 Westendorp 


Plaza et ^|2001b[ [Balthasar & Gomory||2008| ) have been re¬ 
ported, but see also jSchroter \\91\) for smaller vertical field 
gradients of about -0.5 G km For a detailed review of ear¬ 
lier results, see jSolanla] ( [2003 j ) . 

The positive vertical magnetic field gradient observed in the 
middle part of the sunspot penumbra is partly caused by the par¬ 
tial cancellation of the Stokes V signal at the unresolved inter¬ 
face between the spines and the partly, oppositely directed field, 
found at the edges of the filaments, near their heads. However, 
this is also an artefact of the highly corrugated iso-r surface in 
the inner penumbra , see [Balthasar & Gomoryj ( [2008[ ) ; [Balthasar 
[et al.[ ( [2013[ ); [Joshi| ( [2014[ ); [Joshi et al.[ ( [2015[ i. The last two in¬ 
vestigations present a particularly thorough analysis of such ob¬ 
servations. 

The main indication of a magnetic canopy in our inversion 
results is the presence of a stronger field at the upper node than 
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Fig. 15. Left panels: 2D scatter plots of B vs vlos foi' all points within the sunspot boundary at the three optical depth nodes. The 
horizontal dashed line in each panel indicates the zero velocity level. Right panels: the same, but only showing pixels from the 
penumbra of the sunspot. The two horizontal solid lines in panels (e) and (f) indicate thresholds of the velocity used to identify, in 
an approximative manner, the heads and tails of the penumbral filaments. The two vertical dashed lines in panel (f) bound the field 
strength range 1000 < 5 < 2500 G, a condition used to identify heads of penumbral filaments in Fig.p^d). 


at the two lower nodes. Such a signature is only found out¬ 
side the sunspot. This agrees with earlier findings of|Giovanelli| 
T980|;[Gio7anelli & Jones| ( [T98^;[Solanki et al.|(119921 1 19941 


999| ); [Adams et al.| ( |1993 ); Balthasar & Gomory ([2008 1, bu t 
partly differs from the results of jWestendorp Plaza et al.| ( |2001b| ), 
who found a canopy-like structure starting from the middle of 
the penumbra and continuing outside the penumbral boundary 
of their spot. A similar canopy-like structure to that found by 
Westendorp Plaza et al. ([2001 b| ) was repo rted more recently by 
Borrero & Ichimoto| ( |2011j ). Tti^ fact that [Mathew et al.| ( |2Q03| ) 
found no canopy structure, even outside the visible boundary 
of their sunspot, is probably due to the much lower formation 
height of the 1.56 yum lines that they used and, as such, need not 
be inconsistent with our r esults. The sunspot canopie s inferred 
from the 1.56 /tm lines by Solanki et al. ( [1992[p^^ were de¬ 
duced using another technique, and not from the explicit height- 
dependent inversion of spectral lines. Although not so clearly 
marked, a canopy structure outsid e the sunspot boundary , simila r 
to that obtained here, was seen by [Balthasar & Gomory [ ( [20081 ). 


The larger azimuthally averaged inclination of the field by 
about 10° (i.e. on average more horizontal field) in the lower 
layers throughout the spot is con sistent with the findings o f 
[Westendorp Plaza et al.[ ( [2001b| ) and [Borrero & Ichimoto[ ( [2Ql l] ), 
who attributed this to the field canopy structure starting in the in- 
ner penumbra. However, it disagrees with the finding by [Mathew 
[et al.[ ( [2003[ ) in which the field becomes more horizontal with 
height. This disagreement could be caused by the differen t for¬ 
mation height of the IR lines used by [Mathew et al. ( 2003[ ) or by 
the stronger influence of stray light on their results. 

Although we expect the field to be more vertical with height 
in a simple monolithic flux-tube model of a sunspot, a change of 
10° within the lower photosphere is too large to be explained in 
the context of a homogeneous flux-tube model. Instead, this can 
be understood in terms of small-scale structure. Umbral dots and 
light bridges that contain stronger horizontal fields than in the 
umbral background, are visible only in the deepest layers (see 
for example, Riethmilller et al.|2013[ [Lagg et al.|2014 for um¬ 
bral dots and light bridges, respectively, and references therein). 
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but do not contribute to the averaged inclination at higher lay¬ 
ers. Similarly, in the penumbra, the restriction of the mainly 
horizontal-field penumbral filaments to the low photosphere and 
the wrapping of the mag netic field around them ( [Borrero et al.| 
2008} |Tiwari et al.|2Q13| ) leaves only a relatively vertical field in 


the upper layers of the solar photosphere, as proposed by |SolankT| 
|& Moiitavon| ( |1993| ). 

The general inverse relationship between the magnetic field 
strength, B, and its inclination , 7 , is well known from earlier 
observations with lower spatial resolution ( Solanki et al.||1993[ 


An average downflow of about 350 ms“^ is seen over the 
sunspot penumbra at log(T) = 0, caused by the presence of 
strong downfiows at the outer penumbral boundary. The strong 
downfiows are found to continue outside the outer penumbral 
boundary, even reaching a maximum there in the azimuthal av- 


erages (see Fig. [6 

). This is in agreement with the results of, for 

example, [Sheeley 

( I972[);[Dere et al. (1990); Boerner & Kneer 

dl992 

) and Solan 

d et al. ( 1994|), that a part of the mass trans- 


Stanchfield et al.|19W}|Westendorp Plaza et al.|2001b[ [Mathew 


et al.|2004[ ). We revisit this relationship with our high resolution 
data set and find a stronger anti-correlation between these two 
parameters than has been obtain ed by earlier researchers (e.g. 
Westendorp Plaza et al.| ( |200Tb| )). We also find the presence of 


two populations, clearly separated at log(T) = 0: one of umbral 
and of penumbral spine pixels, the other of penumbral filament 
pixels. At higher layers the two populations merge, possibly be¬ 
cause the spines expand and their field wraps around the fila¬ 
ments. 

Another novel feature that we found is a positive correla¬ 
tion between the field strength and the inclination for y > 90°, 
most particularly at log(T) = 0. We attribute this to the tails of 
the penumbral filament s, which contain op posite polarity and a 
stronger magnetic field ( [Tiwari et al.|2013| ). 

It has been known for decades that the d arker umbral re¬ 
gions contain more intense magnetic fields (jKopp & Rabin 


ported by the Evershed flow continues beyond the sunspot’s vis¬ 
ible boundary. 

We can see in Fig.j^that the azimuthally averaged velocity in 
the two top nodes is not zero in the umbra, but instead shows an 
average downflow that increases with height. Since the density 
decreases rapidly with height, any downward mass flow would 
result in such a height dependence of the velocity. This downflow 
could be a result of the inverse Evershed flow ( |Maltby|1975| ). 

The clear pattern of hotter upflows than downflows in the 
sunspot penumbra (by an average difference of about 400 K 
between heads and tails of penumbral filaments) is consistent 
with the transport of heat by convective motions in the sunspot 
penumbra, as has variously been proposed to explain the ob- 


served penumbral brightness 1 

Danielson 19611 [Chitre 

1963| 

Meyer et al. 

1974 

Weiss 

1991 

Jahn & Schmidt[1994[[Sc 

liissler 

& Knolker[ 

2001 

Weiss 

12002 

[Weiss et al.[[20041 [Thomas & 

Weiss 2008 

Rempel et al. 2009b). At the same time, 

, how- 


1 992} [Martinez Fillet & Vazquez|1993[[Stanchfield et al.|1997[ 


Westendorp Plaza et al.[[2001bt [Solanki|[2003[ ). We confirm the 


general relationship at all heights. In qualitative agreement with 
earlier researchers, we find that the darker, and therefore cooler, 
regions contain stronger fields in the umbra at all heights, with 
a dependence roughly in the shape of a raised elephant’s trunk 
(see Fig.[T3]). This is particularly visible in the upper two node 


ever, the downflow speed is correlated with field strength, and 
strong downflows are associated with stronger fields than up¬ 
flows (by an average difference of about 300 G between heads 
and tails of filaments). This result is consistent with the re¬ 
quirements of the siphon flow model ( [Meyer & Schmidt[[I968[ 
Montesinos & Thomas[I997 ). Therefore, the up- and downflows 


in the penumbra exhibit properties that make them partly con¬ 
sistent with both main rivalling theories of the Evershed flow 


positions, indicating that the temperature is nearly constant for a 

et al.[( 1998a|b); Borrero et al. 

([2005 [), and field-free gap mod- 

range of the strongest umbral fields. 

els: [Spruit & Scharmer[([2006|); 

Scharmer & Spruit[p006|)). One 


In the penumbra, however, this relation is not universally 
valid. The heads of penumbral filaments contain strong fields 
( Tiwari et al.[[20I3'] ), yet they are bright, giving rise to a posi¬ 


tive correlation between the temperature and magnetic field in 
the sunspot penumbra. In the early literature, this positive corre¬ 
lation between temperature and field was mis-interpreted as an 
indication that the spines were the brighter par ts of the sunspot 
penumb r a (e.g. Westendorp Plaza et al.[([200Ib[ ), see a review by 


[Solanki[ ( [2Q03[ ), and [Tiwari et al.[ ( [2QI3[ ) for detailed clarifica¬ 
tion). 

The average LOS velocity shows dominant upflows (peak at 

-300 ms“^) in the inner penumbra and downflows (peak at 

1300 ms“^) in the outer penumbra, in qualitative agreement 


with earlier observations (e.g. 

Rimmele|1995bl [Schlichenmaier 

& Sc 

imidt[2000[ [Westendorp 

Plaza et al.[2001a[[Tritschler et al. 

2004 

[Sanchez Almeida et al.[2007t[Ichimoto et al.[2007at [Franz 

& Sc 

ilichenmaier[2009[). This observation supports the scenario 


caveat is that, to test the siphon flow, we need to know v and B at 
a given geometrical height, which is not fulfilled (to an unknown 
extent) by the present data. 

The average azimuthal twist within the sunspot penumbra 
for the deepest layers was found to be -3.8°, which is in agree¬ 
ment with the twist found by [Tiwari et al.[ ( [2009| ) for this spot 
(as observed one day earlier), and estimated by computing the 
spatially averaged signed shear angle, which returns the twist of 
a spot, irrespective of its shape and force-free nature of the field. 
The twist in the spot increases with the radius as well as with 
height. The increase in the twist with height can be understood as 
a result of the expansion of a twisted flux tube that leads to an in¬ 
crease in the azimuthal component of the magnetic field, in turn 
resulting in an increase in the measured twist from tan~^( ^(/,/^^) 
according to |Ven katakrishnan & Tiwari] ( [2009| ); (see, e.g. [Parker 


I974[ [I975[ \\919) . The increase in the twist with radius might 


of the Evershed flow rising in the inner penumbra and sinking in 
the outer penumbra. This average flow pattern is the sum of the 
flows in individual penumbral filaments. The upflows are con¬ 
centrated in the heads, the do wnflows in the tails of individual 
filaments ( [Tiwari et al.|20I3| ). Since the heads are closer to the 
umbra and the tails closer to the outer boundary of the spot, az¬ 
imuthal averaging produces upflows in the inner and downflows 
in the outer penumbra. The upflows a long the axis of th e fila¬ 
ments and downflows along their sides ( [Tiwari et~ar 20I3[ ) partly 


average out during the azimuthal averaging and contribute less 
to the global radial trend. 


be a result of the Coriolis force acting on the outflowing mate¬ 
rial of the Evershed flow. This effect becomes stronger at larger 
radii in the penumbra be cause of the weaker average field (see, 
for example, [Peter| 199^ . A caveat worth mentioning on the de¬ 
termined twist is given by the results of Biihler ( 20I3[ ), who find 
systematic, unexplained positive twists in all features of both the 
magnetic polarities. 

A crude estimate of the mass flux over the full sunspot sug¬ 
gests that the downflowing mass is about 2.5 times larger than 
the upflowing mass. This result is in agreement with that ob- 
tai ned by [Westendorp Plaza et al.[ \\991) for a full sunspot and 
by [Tiwari et al. ' p0I3| ) for a standard penumbral filament. The 
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Fig. 16. Scatter plots of different physical parameters of the standard penumbral filament studied by 


Tiwari et al. 


(2013). The selected 


og(T) = 0. The four 


area also contains pixels belonging to the surroundings of the filament (spines). All plots are for the height 
panels show scatter plots between: (a) the field strength B and inclination y, (b) B and temperature T, (c) the line-of-sight velocity 
vlos and T, and (d) vlos and B. The vertical line in (a) indicates y = 90°; horizontal lines in panels (c) and (d) represent vlos = 0- 
Different parts of the filaments are identified in each plot. 


latter authors suggest that this excess downflow can be attributed 
to the corrugation of the surfaces o f constant optical depth w ithin 
penumbral filaments. In contrast, Puschmann et al. ( |2010| ) find 
the upfiowing mass flux to be five times larger than the down- 
fio wing mass flu x in pa rt of a sunspot penumbra. As described 
in [Tiwari et al.| ( [2013 ), this inconsistency could be because 
[Puschmann et al.||2010| ) evaluate mass flux on the inner part of 
penumbra, where upflows dominate due to the dominant pres¬ 
ence of heads of penumbral filaments, and also partially because 
they take the diskward part of the penumbra, which is blue- 
shifted owing to the Evershed fl ow. Another factor that c ould 
contribute to the difference is that [Puschmann et SI| ( |2010j ) esti¬ 


mate the mass flux on a surface of constant geometrical height, 
which they deduce based on the requirement of divB = 0, 
whereas the ma ss flux was calculated on a sur face of constant 
optical depth by [Westendorp Plaza et al.[ ( [l997[ ) and by us. 


5.2. Constituents of the penumbra 


The filamentary nature of the sunspot penumbra has been inves- 


tigated for more than 50 years (e.g. Danielsoi 

1 ( 1964) 

); Moore 

(1981b|a|); [Title et al. ( 

1993 

); kites et al.[ ( 

19^ ] 

R^immele 

( 1995aj);[Langhans et al. ( 

2005 

); Ichimoto et al. i 

P007b): 

; Borrero 

[et al.[ ([2008[); [Borrero & Ichimoto[ ([2011J; [Joshi et a , 

. (|2011|l; 
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Scharmer et al.|(201 1 ); Scharmer & Henriques| ( [2012| ); Scharmer 


et al.| ( |2013|);|Tiwari et al.| ( |2013| )). One proposal, given by |Lites 


et al.j ( |199^ , divides the sunspot penumbra into two compo 


nents: spines, the more vertical and stronger fields; and inter¬ 
spines: more horizontal and weaker fields. Our results clearly 
show (see Fig. that the spines have properties similar to the 
umbra, but with increasing field inclination and decreasing field 
strength radially outwards. It is worth mentioning that these lo¬ 
cations (umbra and spines) more closely satisfy force-fr ee con- 
ditions on the photosphere than other parts of sunspots (|Tiwari| 

|20T^ . 

The remainder of the penum bra resembles the st andard 
penumbral filament presented by [Tiwari et al. ( 2013 ), who 
show that the penumbral filaments display similar properties to 
stretched granules. They find that penumbral filaments contain 
strong upfiows at their heads and downfiows at their tails. The 
upflows continue along the central axis of filaments to more than 
half of their lengths, surrounded by weak downfiows on both 
the sides of fila ments, signatures of which were already rep orted 
by, for example, jJoshi et al.|(2Qllj ); Scharmer et al.| ( pQTT] ), and 
[Scharmer & Henriquesj ( |2()12| ). The locations of these lateral 
downflows near heads of filaments were found to contain a field 
of opposite polarity to that of umbra, spines, and heads of fila¬ 
ments in one third of the penumbral filaments that [Tiwari et a 


([201 3 [) studied, in agreement w ith the re sults of [Rempel 


fRuiz Cobo & Asensio Ramos [ ( [201 3 [ ) and [S charmer et al. 


( 2012 ); 


(2013 


Additional features, often counted as constituents of the penum¬ 
bra, such as penumbral grains, turned out to be heads of penum¬ 
bral filaments. 

In Fig. we display scatter plots of a number of physical 
quantities at log(T) = 0 from the standard penumbral filament 
created by [Tiwari et al.[ ( [2013| ). By choosing a fixed width for 
all filaments contributing to the standard filament, the latter also 
includes parts of neighbouring spines at its sides and around its 
head. The contributions of different parts of the filament to the 
clouds of data points are marked by arrows. 

A comparison of the scatter plots in Figs.[^a), (b), (c), and 
(d) with Figs[TTJf),[^f),[^f), andp~5]^f), respectively, reveals 
strong similarities, supporting our speculation that penumbrae 
can be described as a combination of spines and filaments. The 
main difference is the lower scatter in Fig. which is due to 
the averaging carried out to produce the standard filament. The 
smoothness, e.g. of spines, and the fact that many points de¬ 
scribe curves in scatter plots of the standard filament, is a result 
of averaging and bicubic s pline interpolation a long the length of 
filaments (for details, see, [Tiwari et al.[[2013| ). In addition, the 
surroundings of filaments in the outer penumbra contain fewer 
spines. Instead filaments lie directly next to each other, thus di¬ 
luting the values for spines and making them smoother in the 
scatter plots. 

6. Conclusions 

Using the SPINOR inversion code, we performed a spatially 
coupled, depth-dependent inversion of a sunspot, observed al¬ 
most at solar disk centre by the Hinode (SOT/SP). We investi¬ 
gated the spot’s thermal, velocity and magnetic structure. The 
average vertical field gradient in the sunspot umbra near the 
deepest node is -1.4 G km“^ which is in agreement with that 
found in the recent MHD simulations. The azimuthally averaged 
magnetic field and inclination show a global general trend as 
identified in earlier work: the field strength decreases with the 
radial distance from spot centre and with height, and the field 
inclination increases with radius, but decreases with height. A 


positive vertical gradient of the magnetic field strength in the 
middle of the penumbra is found that could be either the result of 
a strongly corrugated optical depth surface in the penumbra, the 
signal cancellation of unresolved fin e structure in the deep pho¬ 
tosphere, or a combination of both ( Joshi et ^[2015| ). We also 
find a magnetic canopy structure outside the visible boundary of 
the penumbra. The temperature displays a slightly stronger ver¬ 
tical gradient in the quiet Sun than in the penumbra, but this is 
much stronger than in the umbra. The flatter T(r) profile of the 
umbra is in agreement with earlier investigations (see, for exam¬ 
ple, [SoIanHj5o^ and references therein). 

The strongest fields are found in the darkest regions of the 
umbra, but also in the nearly vertical, opposite polarity flux con¬ 
centrations in the tails of penumbral filaments. Horizontal fields 
in the penumbra show an average field strength of 1 kG, in agree¬ 
ment with the st rength obtained in th e bulk of a standard penum¬ 
bral filament by [Tiwari et al.[ ( [2QT^ . 

The well-known general trend for the magnetic field to be 
strongest in the darkest regions is found in the umbra and in 
spines. However, in penumbral filaments, in particular at their 
heads and in their tails, the opposite trend is found. The field 
azimuth shows an average negative twist of < 5°. The twist in the 
penumbra increases both with radial distance from spot center 
and with height. 

Azimuthally averaged LOS velocities at the deepest layer 
display an upflow (of 300 m s“^) in the inner penumbra, where 
the number of filament heads harbouring strong upflows is 
largest. In the outer penumbra, where the tails of filaments out¬ 
number the heads, an average downflow of 1300 m s“^ is ob¬ 
served. On average, the sunspot penumbra contains a down¬ 
flow of ~350 ms“^ All over the sunspot penumbra, upflows 
are found to be hotter (by about 400 K) than downflows, thus 
qualitatively supporting the thesis that convection accounts for 
the observed brightness of the penumbra and is responsible for 
driving the Evershed flow. However, we also find that strong 
downflows are, on average, associated with stronger fields (by 
about 300 G) than upflows, as required by the siphon flow model. 
Therefore, these data do not in themselves allow us to distin¬ 
guish between these two mechanisms, partly because the actual 
geometric height corresponding to a given r level is unknown. 

A crude estimate of the mass flux over the full sunspot re¬ 
veals that the total downflow of mass exceeds the total upflow- 
ing mass by a factor of ~2.5. This imbalance is probably a con¬ 
sequence of the optical depth corrugation. Since the upflows are 
associated with hotter and less strongly magnetized gas than the 
downflows, we expect that the upflows have been detected at a 
greater height than the downflows. Thus any upflowing mass that 
turns over and flows down again between these two levels is not 
detected, although the associated downflow is seen, giving rise 
to an excess in the downflowing mass. 

We find that the spines in the penumbra display qualitatively 
similar characteristics to those in the umbra. They are, however, 
warmer and have a weaker and more inclined field, consistent 
with their larger distance from the centre of the spot. 

The scatter plots of the different physical quantities of the 
standard filament, including its surroundings, show a qualita¬ 
tive similarity with the scatter plots of the same quantities for 
the full penumbra. In particular, the standard filament shows all 
the populations of points that are also visible in scatter plots 
of the whole penumbra. This confirms that the spines, together 
with the penumbral filaments, are the fundamental constituents 
of the sunspot penumbra. In particular, this similarity suggests 
that there are no further major components of the penumbra be¬ 
sides these two. 
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